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Elastic energy absorption measurements versus temperature on semiconducting, semi-insulating 
(SI) and Fe-doped InP are reported. A thermally activated relaxation process is found only in the SI 
state, which is identified with the hopping of H atoms trapped at In vacancies. It is proposed that 
the presence of In vacancies in InP prepared by the liquid encapsulated Czochralski method is due 
to the lowering of their energy by the saturation of the P dangling bonds with H atoms dissolved 
from the capping liquid containing H2O. The conversion of iron- free InP to the SI state following 
high temperature treatments would be due to H loss with the transformation of the H-saturated In 
vacancies, Vin-H4 donors, into neutral and acceptor Vin-H„ complexes with n < 4. Such complexes 
would produce the observed anelastic relaxation process and may also act as deep acceptors which 
neutralize unwanted donor impurities. 
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I. INTRODUCTION 

Semi-insulating (SI) InP is usually obtained by adding 
Fe before growing the crystal by the liquid encapsulated 
Czochralski (LEG) method. SI InP:Fe offers remarkable 
perspectives as a substrate in optical fiber communica- 
tion systems, and in high frequency electronic devices. 
Essential prerequisites are both high electrical resistiv- 
ity and carrier mobility, in addition to a high material 
homogeneity. Generally, the semiconducting as-grown 
InP is quite inhomogeneous, as it exhibits growth stri- 
ations, dislocations decorated by microdefects, micro- 
precipitates embedded in the crystal lattice and some- 
times inclusions.tl In particular, since the Fe concentra- 
tions necessary to reach the semi-insulating state are 
close the Fe solubility limit, domains of Fe precipitates 
are inevitably formed. Post-growth thermal treatments 
of the ingots at about 900 °G produce beneficial effects 
concerning defect annealing besides a resistivity increase, 
but cause inhomogeneities across the cwstal diametercl 
and do not improve the Fe homogeneity.cl Therefore, the 
thermal treatments have to be carried out on the wafers 
cut from the as-grown ingots; in this way the wafer ho- 
mogeneity results improvedu and the striations and other 
types of lattice imperfections may be annealed. 

It is diffusely believed that in the SI state of InP:Fe, the 
electrons supplied to the crystal by the unwanted residual 
impurities (shallow donors) are trapped by Fe which acts 
as a deep acceptor and provides compensation. However, 
the SI state is reached only if the Fe concentration (typ- 
ically in the range of 10^^ cm^'^) is remarkably higher 
than the residual donor impurity content; as a confir- 
mation, it was reported that the impurity content hi SI 
InP:Fe is lower than the free carrier concentrationjj in- 
dicating that sources of donors other than the impuri- 
ties considered above must be actb/e; they were identi- 
fied later as Vin(HP)4 complexes Jj'lI and labelled here as 
Vin— H4. Lastly, it should be mentioned that the achieve- 
ment of the SI state by Fe doping may present further 
disadvantages due to the nonuniformity of the electrical 
properties along the crystal growth axis caused by Fe 
segregation.! 

The semi-insulating state is also reached by thermally 
treating undoped InP for several hours or days at about 
950 °G under phosphjomus atmospheres ranging from 
tenths to tens of barsj3"llil after which the resistivity can 
be increased above lO"" Q cm. However, for the conversion 
to occur it is necessary that the carrier concentrations be 
below the limit of 4x10^^ cm~^ (Ref. |lO|); thermal treat- 
ments of materials having higher residual carrier concen- 
trations due to higher impurity contents in the starting 
material usually display only a slight reduction in the 
concentration of free carriers after annealingji3 and the 
same effect is also produced in SI InP:FeJa 

Several mechanisms have been proposed for the SI con- 
version, and many of them hypothesize that the thermal 
treatment may activate some kinds of deep acceptors in 



such a concentration to compensate the shallow donors 
related to impurities. Among the possible deep accep- 
tors, the Fe atom both in Fe-doped and in undoped InP 
(therefore present as residual impurity) has been con- 
sidered. However, also the hydrogen-related defects are 
receiving great attention, as H is always presentllj and 
may form complexes with native defects and impurities. 

Generally, the studies of the properties of SI InP 
are conducted on the material doped with Fe, whilst 
a systematic analysis of the comparative properties of 
SI InP:Fe and undoped InP is missing. The purpose 
of the present work is to carry out, through measure- 
ments of the complex modulus (elastic energy loss and 
dynamic modulus), a comparative study of the defects in 
InP in the semiconducting and the semi-insulating state 
obtained both by the doping with Fe and by the thermal 
treatment. 



II. EXPERIMENT 

The thermal treatments (TT) finalized to the conver- 
sion from the semiconducting to the SI state were carried 
out on wafers at about 950 °C under a phosphorous vapor 
pressure of 0.2 bar for times of the order of hours. 

The crystals were grown by the LEG method at 
Venezia Tecnologie, and the samples were three rectan- 
gular bars (about 43 x 5 x 0.4 mm^) of normal purity and 
high purity. They were characterized as follows: sample 
1 was an undoped highly pure sample subjected to the 
TT for conversion and presented the SI state. Sample 2 
was of normal purity and doped with 10^^ cm~^ Fe impu- 
rities; therefore it was therefore semi-insulating without 
need of TT. Sample 3 was highly pure but not subjected 
to the thermal treatment for the conversion, and hence 
non semi-insulating. 

The elastic energy loss and resonant frequency mea- 
surements were carried out between 1.9 and 650 K. The 
samples were horizontally suspended on two nodal lines 
by thin wires and three flexural vibration modes were 
electrostatically excited in the frequency range between 
1 and 16 kHz. The vibration amplitude was detected 
through the variation of the sample-electrode capacitance 
by a frequency modulation technique, the elastic energy 
absorption was measured from the rate of the free decay 
of the vibration. 

The elastic energy loss coefficient, or reciprocal of the 
mechanical Q of the sample, is related to th&|imaginary 
part of the complex Young's modulus E* asEJ 

Q-^ - l^r. ('-§-) ^E^iAXf-^, (1) 
\ E J k^T 1 + (tjr) 

where the last expression is the contribution of a molar 
concentration c of defects diffusing or reorienting with 
a characteristic time r between states or sites which are 
made inequivalent by the vibration stress cr; AA is the ap- 
propriate component of the change of the elastic dipole 
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tensor after a defect reorientation, and the elastic dipole 
tensor A is defined as the strain caused by a unit concen- 
tration of the defects; vq is the cell volume and uj — 27r/ 
is the angular vibration frequency. The above contribu- 
tion to the absorption is peaked at the temperature at 
which ujT = 1 ; since r is temperature dependent and usu- 
ally follows an Arrhenius law, r (T) = tq exp (iJ/fceT), 
the peak shifts to higher temperature if measured at 
higher u;. The present samples had the (110) axis par- 
aUeJ to the length, so that the appropriate modulus 
Si2 + ^544)- For trigonal defects 

A2) while for (110) or- 
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thorhombic defects 
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S12) + 1^544 — 

The Young's 



(Ai + A2-2A3)" + (Ai-A2) 

modulus is related to the frequency of the first flexural 
mode by = p (0.975 / P/t)^, where p, I and t are the 
sample density, length and thickness; for our samples it 
is E[iio] = 1.40 X 10" Pa. 



III. RESULTS 

Figure 1 shows the elastic energy dissipation curves 
from 1.9 to 370 K of the undoped SI sample 1 and sample 
2. An intense peak (peak A) is present at about 290 K 
in both samples, whereas the dissipation curves in the 
low temperature side indicate that two additional peaks 
(around 130 and 220 K) may be present but confused 
with the background; however, due to their small height, 
they will not be considered here. It is surprising that 
the two samples, having in common only the SI state but 
not the purity and the thermal history, present the same 
type of spectrum, and differ only for the height of the 
background dissipation. 

In order to test the thermal stability of the anelastic 
relaxation spectrum, sample 2 was annealed in vacuum 
for about 30 min at temperatures subsequently increased 
from 400 to 660 K at steps of 40 K; after each step the 
relaxation spectrum was measured, and the results are 
reported in Fig. 2. With increasing annealing temper- 
ature the background dissipation progressively decreases 
down to a saturation value which is attained after the 
treatment at 620 K (7'*^). Instead, the height of the peak 
remains practically unaffected, and is also stable versus 
thermal cycling, as the curves are retraced on heating and 
cooling. After this series on annealings, it was checked 
that the sample remained in the semi-insulating state. 

The energy-loss curves, measured at three vibration 
modes (between 1 and 16 kHz) showed that the peak A 
is thermally activated, since it shifts towards higher tem- 
peratures with increasing frequency. Figure 3 presents 
the experimental data of sample 2 after the 7"^ vacuum 
annealing at 620 K at 1.2 and 15.6 kHz, after subtrac- 
tion of the background; the data are satisfactorily fitted 
by Debye curves (continuous lines) with a single relax- 
ation time following the classical Arrhenius law. The 



activation energy of the relaxation process and the pre- 
exponential factor of the relaxation rate obtained from 
the fit are E = 0.67 eV and t^^ = 3 x lO^-* s"^ 

The spectrum of the semiconducting sample 3 is re- 
ported in Fig. 4 together with those of samples 1 and 
2 (continuous lines) for comparison. Sample 3 displays 
(curve 1), superimposed to a high background, a small 
peak at 220 K and a well developed peak at 400 K which 
decreases during its measurement and disappears after a 
vacuum annealing of about 30 min at 450 K (curve 2). 
It is remarkable that the peak A, found in the SI sam- 
ples, is absent here, or at least non emerging from the 
background. The lack of peak A in the non-SI samples 
was confirmed also in undoped InP samples cut from the 
same ingot and from a different ingot (results not shown) . 

In order to study the stability of the anelastic spectrum 
in the non-SI sample 3, a series of vacuum annealings of 
30 min each at increasing temperatures was conducted, 
as for sample 2. The results are displayed in Fig. 5 and 
show that: i) the dissipation around 400 K is not a sta- 
ble peak as a function of temperature; indeed, it starts 
decreasing on heating above 350 K, even keeping the sam- 
ple temperature constant (result not reported here) and 
practically disappears after a vacuum annealing at 500 K; 
ii) the small peak around 220 K increases after the first 
vacuum annealing (450 K), but is depressed after the 
subsequent treatments at higher temperatures; Hi) the 
overall dissipation steadily decreases down to saturation 
values with the subsequent isochronal annealings at pro- 
gressively higher temperatures. 

Figure 6 compares the dissipation curves of SI InP:Fe 
(sample 2) and of non-SI InP (sample 3) after the last an- 
nealing treatments at 650 K; only the peak at 300 K in 
the SI sample persists after the high temperature treat- 
ments. 

In order to test whether peak A is suppressed by heat- 
ing above 660 K, the isochronal annealings on sample 2 
were continued up to 847 K, and Fig. 7 shows that the 
species giving rise to peak A is stable up to that temper- 
ature. At this temperature a few droplets appeared on 
the sample surface, indicating surface degradaticua-due to 
P loss and consequent formation of In clusters .Il3 How- 
ever, the sample bulk was very little affected, as deduced 
by the small relative variation of the resonance frequency 
(< 3.5 X 10"'^); the vibration frequency is directly related 
to the elastic modulus, which is very sensitive to the for- 
mation of second phases. In conclusion, peak A is stable 
up to the sample decomposition in vacuum. 

To better interpret the annealing behavior observed, 
two samples cut from the same wafers of sample 2 (Fe- 
doped) and sample 3 (non-SI) were subjected to vacuum 
effusion experiments with a rate of 3-4 K/min, monitored 
with a residual gas analyzer. It was found that H out- 
diffusion starts above 560 °C, i.e. at about the same 
temperature at which the sample decomposition starts, 
as indicated by the appearance of some peaks in the ef- 
fusion spectrum corresponding to the masses of PH„ and 
P2H„, besides H2. 
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IV. DISCUSSION a nearest neighbor bonded to P. Peak 2 is of uncertain 

origin, but fflay be due to a cluster similar to that of peak 
A. The vacancy-H complexes 3 with MnO 



Optical absorption experiments on InP:Fe have re- 
vealed several local vibrational modes (LVML regard- 
less of whether—the material was aSigmwnJjO plasma 
hydrogenatedjQO or ion implantediEfTcEl these lines oc- 
cur in the range included between 2200 and 2300 cm^^ 
of the spectral region. Howevet- if the same mate- 
rial is implanted with deuteronsc2l instead of protons, 
the corresponding LVM are shifted to lower energies 
(~ 1600 cm~^), and the ratio between the energies is 
close to (ttid/toh)"'^^^- This observation provided evi- 
dence that the peaks are due to hydrogen related defects 
cause, apjd—thBy have been attributed to PH stretching 
modes.BEiy-E2r 

Four of the reported lines are particularly important, 
as they have been reproducibly observed, and there seems 
to be an agreement on their attribution; such peaks occur 
at 2202 cm"^ (peak 1), 2273 cm'^ (peak 2), 2286 cm^^ 
(peak 3), and 2316 cm~^ (peak 4). The main line is peak 
4 (Ref. but peak 1 is oLcomparable intensity if the 
samples are ion implantedpj instead, peaks 2 and 3 are 
rather weak. 

Peak 1 is originated by a H complex with trigonal 
symmetry^ jsszith a P-H bond oriented along the (111) 
dircction;LlO the low energy of its stretching mode is 
indicative of a long bond length: therefore it has been 
attributed to an In vacancy having one the P dangling 
bonds saturated by a H atom occupying a bond centre 
(BC) site between Vin and P; this complex will be la- 
belled as Vin-H. 

Peak 4 appears in the high energy region of the PH 
stretching mode, indicating a shorter bond length. Op- 
tical absorbance measurements in undoped InP samples 
cocharged with H and D have convincingly demonstrated 
that the complex giving rise to this peak is constituted 
by more than one H atom.Q In addition, the stress split- 
ting of this line indicates a cubic symmetry for the cor- 
responding cluster .Q The various authors agree that the 
complex consists of an In vacancy decorated by four H 
atoms in BC occupancy, Vin— II4; in this configuration, 
the repulsion between the four H atoms shortens the P- 
H bond lengths and gives rise the observed high energy 
LVM. The four H atoms, each bonded to one of the four 
P atoms around Vin, contribute four electrons to the 
trivalent vacancy and leave one electron available .Q The 
Vin— II4 complex is therefore a shallow donor, and could 
be responsible for the discrepancy between the donor con- 
tent estimated from the concentrations of the impurities 
and f]pOm the Fe+ content that is necessary to passivate 
them.cl 

Peak 3 was studied by local vibrational mode ex- 
periments i-#J^d was attributed to the Zn residual 



impurities ;EI[£!J the complex is believed to be formed by 
the Zn atom occupying an In site having one H atom as 



B. Peak A 

In InP, the omnipresence of hydrogenEl (certainly due 
to the LEG preparation method) , the fact that LVM lines 
due to isolated H hajjejipt been reported, the ability of H 
to form complexesBLlllj^EJ and its high mobility, strongly 
indicate that the atomic defect giving rise to peak A may 
be a H-related complex. Undoubtedly, the presence of the 
thermally activated process with maximum near 300 K 
proves that the InP samples in the SI state contain a 
species which is highly mobile at room temperature (~ 
lO'^ jumps per second) and that this species is rather 
abundant, considering the height of the corresponding 
peak. In addition, the value of the pre-exponential factor 
is typical of point defect relaxations and indicates that 
this imperfection involves an atomic complex. 

It is remarkable that the process is characterized by a 
single relaxation time: in fact, the absence of peak broad- 
ening indicates that the mobile species relaxes with only 
one type of elementary jump, and that there is no in- 
teraction among the relaxing units. Single-time Debye 
processes are not common, andpXine was also discovered 
in Si:B charged with hydrogen.L3 In that system the H 
atom hops among the four equivalent BC sites coordi- 
nated with substitutional B, overcoming an energy bar- 
rier {E = 0.22 eV) which is lower than that presently re- 
ported. The absence of broadening of the peak due to the 
reorientation of the H-B pair led to the conclusion that 
the interaction between H-related defects may be negligi- 
ble at dopant concentrations as high as 10^^ at/cm^. The 
relaxational dynamics of H within the H-B pairs in Si has 
also been measured by the decay of the stress- induced IR 
dichroismc3 for relaxation times between 10'^ and 10^ ^ 
By joining the data for the two types of experiments, C3l£| 
relaxation rates spanning 11 orders of magnitude were 
obtained in the temperature range 60-120 K; a weak devi- 
ation from the Arrhcnius dependence was observed at low 
temperature, which was interpolatedcj with a polaron- 
like model. 

Similarly to the case of B-H in Si, it would be natural 
to think of the reorientation of H around Fe as a source of 
anelastic relaxation in the Fe-doped InP sample. Indeed, 
measurements of the local vibrational modes show that in 
InP:Fe, H resides between the Fe atom (which substitutes 
In) and a nearest neighbor P atom, mainly bonded to P. 
However, the fact that both undoped and Fe-doped InP 
display the same peak (see for instance Fig. 1, where the 
small temperature shift between the two peaks is due to 
the different frequencies) demonstrate that Fe does not 
play any role in the mechanism causing peak A. 

We discuss now the possibility that peak A is due to 
Vin— H„ complexes (n = — 4). For n = and n = 4 the 
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complexes have in principle the same cubic symmetry of 
the crystal, so that there would not be any differentiation 
between complexes under stress, unless Jahn- Teller type 
lowering of the defect symmetry occurs. Indeed, it has 
been proposed that in certain ionization states the Vin 
is distortedEJ In this case, different distortions would be 
distinguished by stress and would cause anelastic relax- 
ation, but the reorientation between them would be ex- 
pected to occur with a faster rate than in the cases where 
atomic hopping is ipvolved and produce a peak at much 
lower temperature:L3 the barrier of 0.67 eV seems too 
high for the reorientation of a Jahn- Teller like distortion. 

Instead, the reorientation of the Vin— II„ complexes 
with n = 1,2 and 3 is expected to appear as a slower 
single-time relaxation process. In fact, in all three cases 
the complexes can reorient among equivalent configura- 
tions through a single elementary jump of an H atom 
between the four different P dangling bonds around the 
vacancy. For n — 1 there are four equivalent BC sites 
for H, each coordinated with one of the four P bonds, 
and the defect has trigonal symmetry. For n = 3, one 
of the four BC sites is unoccupied, so that the complex 
reorientation is due to the jump of the H vacancy instead 
of the H atom, and the symmetry is again trigonal. For 
n — 2 the H atoms can occupy six equivalent pairs of BC 
sites along the [110] directions. The corresponding elas- 
tic dipole has the principal axes along two (110) and one 
(100) directions and its symmetry is (110) orthorhombic. 

The three complexes may have different energy barri- 
ers for reorientation, producing distinct relaxation peaks. 
In this case, peak A would be due to the dominant com- 
plex, and the other two complexes would produce peaks 
at different temperatures, e.g. the two small peaks at 130 
and 220 K. It is also possible that all the complexes have 
nearly the same energy barrier for reorientation. In fact, 
the LVM energies of the various complexes differ by only 
5% (Refs. ^,0), and this can be taken as an indication 
that the potential felt by H and therefore the barriers 
between different minima are little affected by the pres- 
ence of neighboring H atoms. In this case all the three 
complexes would contribute to peak A, and indeed the 
present data are compatible with a spread of activation 
energies of the order of 5%. At present it is not possible 
to decide between the two cases, and further measure- 
ments on crystals with different orientations are needed. 

Finally, we note that the concentration of defects pro- 
ducing peak A estimated from the peak intensity is fully 
consistent with the concentration of H and Vin-II„ com- 
plexes which are generally present in InP crystals grown 
by LEC,EI i.e. of the order of 10^^ cm^^. in fact, the 
peak intensity in Eq. (1) is 

Q^l. = ^^AA^ (2) 

and substituting the measured E — 1.40 x 10^^ Pa and as- 
suming thattll 0.1 < A A < 1, the volume concentration of 
defects turns out of the order of c/vq 6 x 10^^^^ cm^'^. 



C. The background dissipation and the peak at 
400 K 

The high background dissipation of both the SI (Fig. 2) 
and the semiconducting (Fig. 5) samples might be due to 
the thermoelastic effect or to the movement of disloca- 
tions. The thermoelastic effect is due to the periodic heat 
transfer througl^-the sample induced by the alternating 
vibration stress;til the corresponding relaxation rate 
in Eq. (1) is directly proportional to the thermal con- 
ductivity, which is generally a slowly varying function 
of temperature, so that a very broad peak results. The 
observed lowering of dissipation after the high tempera- 
ture treatments should be mainly explained in terms of a 
change of the charge carrier concentration and mobility, 
which in turn affects the thermal conductivity. During 
the present experiments, however, no relevant changes of 
the electrical conductivity were observed. Therefore, it 
is unlikely that the main contribution to the background 
dissipation in the as-grown samples was of thermoelastic 
origin. 

Dislocations are formed during the crystal growth and 
subsequent cooling, and their stress-induced movement 
may produce a broad background dissipation. The de- 
crease of such a dissipation after the high temperature 
treatments can be caused by the migration of defects 
onto the diskication lines, resulting in the pinning of the 
dislocations.E^ The migrating defects should be rather 
mobile, since the lowering of dissipation starts already 
a few tens of kelvin above room temperature; they can 
be oxygen or hydrogen atoms. It should be noted that a 
few at ppm impurities are sufficient to saturate the whole 
dislocation net. 

The presence of the peak at 400 K in the semiconduct- 
ing sample 3 and its disappearance during its measure- 
ment, which occurs at a temperature at which there is 
certainly no H loss, indicates that the sample is not in 
thermal equilibrium. The species giving rise to the pro- 
cess is again rather mobile; therefore, the peak is likely 
due to a H- related defect, but it is not possible at present 
to make any attribution. 

D. Interpretation of peak A and proposition of a 
model for the semi-insulating state 

From the above results it appears that peak A is 
present only in the SI state of both Fe-free and Fe-doped 
InP, suggesting that there is a common mechanism of 
conversion to the SI state, which is also responsible for 
the anelastic relaxation process. 

The SI conversion in InPiFe is attributed to the re- 
moval of the shallow donors by the deep acceptor Fe 
atoms. Similarly, during the high temperature treat- 
ments in P atmospheres necessary to reach the SI state in 
undoped InP, the activation of some sort of deep accep- 
tors able to compensate the shallow donors introduced 
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by impurities has been unanimously invoked. However, 
there is no agreement at present on the nature of such 
acceptors and consequently on the SI conversion mech- 
anism. Bliss et alB reported that the 2316 cm~^ line 
attributed to Vin— H4 disappears after the TT finalized 
to the SI conversion, and in addition they observed that 
Fe and Cu were introduced during their TT. Op-this basis 
they concluded, in accordance with Hirt et aL,E3 that the 
high resistivity values are due to the contamination by Fe 
and Cu deep acceptors during the annealing process and 
to the concomitant outdiffusion of the Vin— H4 complexes 
acting as intrinsic donors. We exclude that Fe plays such 
a role in our undoped SI samples, since according to the 
chemical analysis sample 1 contained less than 1 atomic 
part per billion, and the TT for the SI conversion was 
carried out under controlled purity conditions. 

For the model we presently propose we assume that 
hydrogen and vacancies are abundant, even though not 
intentionally introduced. As discussed in the follow- 
ing paragraph, we suppose that the as-grown material 
mainly contains Vin— II4 complexes, whereas lower or- 
der complexes are also present if H is trapped by other 
defects, like Fe. Xte assumption is corroborated by sev- 
eral observationsQ'Ea of the LVM at 2316 cm^^ (Vin— II4) 
in as-gro-pm InP, and of the 2202_cm^^ (Vin— H) line in 
Fe-dopedEa and proton-implantcdEHl laP. We believe that 
the disappearance of the Vin— II4 lincO consequent to the 
TT at 950 °C for the SI conversion is due to the outdif- 
fusion of H, rather than of the whole Vin— II4 complex, 
since we expect that Vin is less mobile than H. As a con- 
sequence of the H loss from the sample, transformation 
of Vin— II4 complexes into complexes with n < 4 takes 
place. Indeed, after thermal annealings above 600 °C 
new lines appear at lower-energyjj attributed to Vin— Hn 

with n < 4. nnn 

Considering that Vin is a triple acceptor ,l30'E2l we 
expect that the Vin— II3 complex is neutral, as the H 
atoms supply the three electrons. Similarly, Vin— II2 and 
Vin— H are expected to be acceptor complexes, and we 
hypothesize that they are deep acceptors. The dissocia- 
tion of the Vin— II4 complexes during high temperature 
annealing reduces the number of shallow donors, and the 
newly formed Vin— II„ (n = — 2) deep acceptors remove 
the free carriers introduced by impurities and render the 
material semi-insulating. 

The proposed model is congruent with the observation 
that peak A (attributed to Vin— II„ with n = 1 — 3) in 
undoped InP appears only after the TT at 950 °C, which 
transforms the complexes with n = 4 into lower order 
complexes; therefore, the existence of the peak should be 
closely related to the SI state. Accordingly, we expect 
that peak A is absent in as-grown undoped samples, as 
experimentally observed (Fig. 4). The occurrence of the 
peak also in the Fe-doped sample, which did not undergo 
the TT, can be explained in terms of partial trapping 
of H in Fe-H pairs, with consequent transformation of 
Vin— Il4 centres into lower order complexes. 



E. Formation of V-H„ complexes 

It remains to be explained why InP prepared by the 
LEC method contains a substantial concentration of V- 
H„ complexes, as indicated by the corresponding LVM. 
In fact, the equilibrium concentration of vacancies in 
crystals are very small even near the melting tempera- 
ture, and for the case of InP it seems natural that va- 
cancies of the more volatile P are predominant over the 
Vin. However, the formation energy of Vin can be low- 
ered if the P dangling bonds are saturated by H atoms, 
and could even become negative; then, in the presence 
of H the formation of In vacancies saturated by H would 
be possible or even energetically favorable. During the 
LEC crystal growth, hydrogen woulcL-be supplied by the 
capping liquid, which contains H2 0;E1I H could dissolve 
in the liquid InP forming P-H bonds and then be incor- 
porated in the crystal as Vin-H4. 

That the P-H bonds are much more stable than the 
In-H bonds is demonstrated by the lack of In-H LVM 
in InP even after the irradiation with H, which produces 
both P and In vacancies and self-interstitials. Within the 
proposed picture, the high stabihty of the P-H bond is 
consistent with the stability of peak A against vacuum 
annealings up to the sample decomposition (Fig. 4) and 
with the relatively high barrier for the hopping of H be- 
tween different P bonds in the Vin (0.67 eV, compared to 
0.22 eV for the reorientation of the B-H complex in Si). 

The formation of superabundant vacancies iSpsJready 
observed in metal hydrides at high temperatures J23 when 
the binding energy of H to the metal vacancies can sub- 
stantially lower the V-H„ formation energy. By exposing 
the metal to a high H pressure at high temperature, it 
is possible to let very high concentrations of H-saturated 
vacancies to migrate into the bulk (even tens of percent). 

V. CONCLUSION 

Anelastic spectroscopy experiments on InP crystals 
grown by the liquid encapsulated Czochralski method 
have revealed the presence of defects whose reorientation 
rate is thermally activated over a barrier of 0.67 eV. The 
corresponding anelastic relaxation process is observed 
only in the semi-insulating InP samples, irrespective of 
whether the semi-insulating state is obtained by dop- 
ing with Fe or by annealing near 900 °C. It is proposed 
that these defect are In vacancies with the P dangling 
bonds partially saturated by H, whereas vacancies com- 
pletely filled with four H atoms do not cause anelastic 
relaxation. A possible mechanism for the formation of 
the Vin-H4 complexes and for the reversion to the semi- 
insulating state after annealing near 900 °C is proposed, 
based on the stability of the P-H bonds. The formation 
of such bonds with H dissolved from the capping liq- 
uid during the crystal growth would make the formation 
of Vin-H4 complexes favorable.. Such complexes would 
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be donors which contribute to make the crystal a semi- 
conductor. The subsequent treatment at 900 °C would 
partially outgas H from the sample and remove it from 
the Vi„-H4 complexes, transforming them into neutral or 
deep acceptor defects. 
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VI. CAPTIONS TO FIGURES 

Fig. 1 Elastic energy loss versus temperature of semi- 
insulating InP. Sample 1: highly pure and subjected to 
the thermal treatment for conversion to the SI state; sam- 
ple 2: doped with 10^^ cm'^ Fe. 

Fig. 2 Peak A of sample 2 during a series vacuum an- 
nealings for about 30 min at temperatures subsequently 
increased from 400 and 660 K at steps of 40 K. 

Fig. 3 Peak A in sample 2 after the 7**^ vacuum an- 
nealing at 620 K, after subtraction of the background. 
The continuous lines are single time Debye curves with 



E = 0.67 eV and Tn 



3 X lO^"* s-i. 



Fig. 4 Elastic energy absorption of semiconducting 
sample 3 (/ = 2 kHz) during two subsequent measure- 
ments (curves 1 and 2). For comparison, the dissipation 
curves of the SI samples 1 (2 kHz) and sample 2 (1.2 kHz) 
are also shown. 

Fig. 5 Evolution of the anelastic spectrum of semicon- 
ducting InP (sample 3) measured in vacuum: as prepared 
(curve 1); 2"^^ heating run (curve 2); after subsequent 
annealings in vacuum for 30 min at the temperatures in- 
dicated in the figure. 

Fig. 6 Elastic energy loss of SI InP:Fe (sample 2) and 
of non-SI InP (sample 3) after the last vacuum annealings 
at 660 and 650 K respectively 

Fig. 7 Peak A in sample 2 after annealings of 30 min in 
vacuum at the increasing temperatures indicated in the 
legend; these measurements follow those of Fig. 2. 
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Fig. 7 
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